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Abstract 

The relation between the Poisson and Schrodinger equation is obtained through a simple transformation. 
It is pointed out that this analogy between both equations can be only applied for potentials that involve 
a combination of attractive and repulsive delta function potentials. This relationship enables us to use 
\ elementary electrostatic results to find the ground state energy and wave function of the associated quantum 

state problem in one dimension. Particularly, the result shows that it is possible to trap a single electron in 
a one dimensional ionic crystal modeled by repulsive and attractive delta functions. 
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I. INTRODUCTION 



It is very well known that Poisson's equation is among the most important equations in electro- 
magnetic theory. This equation fundamentally links the source of the field, i.e. charge distribution 
p, and the electric potential V . In theory, this equation allow us to find all we need to know about 
an electrostatic field given the charge distribution. The general solution to Poisson's equation has 
the following integral form' 

= [ T^.d^r'. (1) 

47reo J \r — r'\ 

Once we know V the electric field is calculated from the gradient of V , i.e. E = —W. Nev- 
ertheless, the integral given in Eq. ([T]) can only be evaluated analytically for the simplest charge 
configurations. This is why it is very useful to recast the problem using Poisson's equation with 
appropiate boundary conditions. For the restricted case of static charges and fields in vacuum, the 
appropiate boundary conditions for the electric potential are^ 
/■ 

V is always a continuous function across any boundary, and 

iA{dV/dn) = — cr/eo the normal derivative of the potential is discontinuous at the boundary, 
where a represents the surface charge distribution at the boundary. 

(2) 

On the other hand, the standard boundary conditions for the wave function ^ in one dimensional 
quantum mechanics are :- 

^b is always continuous, and 

(3) 

dip/dz is continuous except at points where the potential is infinite. 

In particular, if there is a delta function potential at zq, then the change of slope of the wave 
function, A{d'ip/dz), at the point zq is proportional to the amplitude of the delta function as well 
as the wave function (/'(-Zq) at this point.- This is the type of boundary conditions we expect for the 
electrostatic potential. As a consequence, the electrostatic potential and the one dimensional wave 
function share the same boundary conditions for a delta function potential. 
The one dimensional delta-function potential has proved to be a useful model to represent diatomic 
ion and molecular orbitals,^ two electron systems,- electron-dipole interaction,^ solids,-^ the many 
body problem,- and scattering.— It is well known that the single attractive delta function potential, 
i.e. U{z) = —aS{z), possesses one exponentially localized bound state for all values of the 
strength of the well a > 0. Its existence and stability has been tested for the effects of different 
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boundary conditions^ and symmetry-breaking perturbations.— 

In this article, a mapping between a three dimensional electrostatic problem described by Poisson's 
equation and the one dimensional Schrondinger equation is given through a simple transformation. 
This relation will allow us to find the wave function and energy of the ground state of the quantum 
system by using electrostatic results. 

The article is organized as follows. First we will start by transforming the Poisson equation into 
a Schrodinger-like equation, and discuss when both equations are equivalent. Then we will apply 
electrostatic results to solve the Schrodinger equation for the case of a potential function given by 
a combination of attractive and repulsive delta functions. The conclusions are summarized in the 
last section. 



II. DUALITY BETWEEN POISSON AND SCHRODINGER EQUATION 

The Poisson equation in three dimensions is given by 

VV = — ^, (4) 

where V is the electric potential, p is the volume charge density and cq is the electric permitivity 
of free space.-^ 

Making the following transformation V{r) = Vo\n{^{r)/A), where Vq and A are constants to 
ensure dimensional consistency, we obtain 

where we have used the fact that E = —VV = —Vq^. Multiplying Eq. ^ by Voeo\E'/2 we end 
up with 

-Yt^^^^-^^Vop^ = -'jE.E^. (6) 

Equation ([61) is a Schrodinger-type equation when the magnitude \E\ is constant. We know from 
basic electrodynamic courses that the only charge distribution to produce uniform electrostatic 
fields are infinite charge sheets. Since the infinite charge sheet has plane symmetry, the electro- 
static potential depends on one variable only, i.e. \1' = ij{z). For this case Eq. © is given by 

Vo^eod^^ 1 €0 2/ 



Making the following substitutions: 

Vo'eoal = ^, -Wop{z)al+'j{El-E'^Jal = U{z) and ieoi?L«o = 1^1 (8) 

where oq is an arbitrary length and Ezoo is the asymptotic magnitude of the electric field, i.e. 
Ezoo = lim^_^oo \Ez\,we obtain the time independent Schrodinger equation 

h"^ cPih 

-—-p^ + U{z)^ = -\E\,p, (9) 

for the case of bound states, i.e. < 0.— This relationship means that the wave function and the 
bound state energy will be given in terms of the electric potential and the energy density by 

iIj{z) = Ae^(")/^«, E = -U, (10) 

where A is the normalization constant and U = eoE'^/2 is the energy density, respectively. Note 
that the wave function has no nodes, therefore only the ground state can be found with this method. 
This is consistent with the electrostatic uniqueness theorem which guarantees that the electrostatic 
potential is uniquely determined for a given charge density distribution and the electrostatic po- 
tential at all boundaries are specified.- 

In order to have physically realizable quantum states the wave function must be normalized . ^^^^ ^ 
In the following section we solve the Schrodinger equation using electrostatic solutions. 

Ill, SOLVING THE SCHRODINGER EQUATION VIA ELECTROSTATIC SOLUTIONS 

We know from basic electrodynamic courses that the only charge distribution to produce a uni- 
form electrostatic field are an odd number of parallel infinite sheets with opposite electric charge 
densities ±a^-^, hence we will consider a volume charge density given by 

N 

p{z) = (T - ^n), for iV = 0, 1, 2, . . . (11) 

n=-N 

where 5{z — Zn) is the Dirac delta function at the position Zn^^ Let us consider the case when 
= 0, for this case the volume charge distribution is of the form 

p{z) = a5{z), (12) 

where we have taken zq = The electric field caused by an infinite sheet of charge density a at 
any space point is^ 



therefore, the electric potential and the energy density are given hy^ 

V{z) = -^\z\, U = ^. (14) 
Using equation (flOl) we obtain the wave function and the energy of the bound state 
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The wave function given by equation (115]) can be normalized only for a > 0, otherwise it blows 
up when z — t- ±oo. Using equation ([8]) we get the potential energy function U{z) = —a5{z) 
for the quantum mechanical case, where a = <jVq/2 is the strength of the potential. Making the 
appropiate substitutions into equation (fT5l) we get the corresponding bound state solution in the 
presence of an atractive delta function potential 

ilj{z) = Ae-^"\'\"^\ E = -^, (16) 

where the normalization constant is A = y/ma/h? Using the wave function given in equation 
(fT6l) it is easy to obtain the expectation values for the potential energy {U)q = —ma^/li? and 
kinetic energy (T)o = mc? j2W', respectively. 

If we consider the general case for > 1, then the electrostatic potential and energy density are 
given by 

N 2 

Interestingly, there is one quantum bound state with exactly the same energy for the charge distri- 
butions given by equations (fTTI) and (fT2l) . where the wave function is given by 

Remarkably, the solution given in equation (fTSi) is the product of growing and decaying exponen- 
tial functions. In order to have a physically realizable quantum bound states we need to normalize 
the wave function given in equation (fTSi) . For simplicity, I am going to consider the case consisting 
of evenly spaced infinite charge sheets, i.e. Zn = na, where a is the distance between adjacent 
sheets. For this particular case the wave function is an even function and in order to normalize ^ 
we need to solve the following integral 

2\A\^ f exp [ ^ ^ (-l)"|2-na| ] = 1. (19) 







n=—N 



Using the fact that 



N 



N 



J2 {-ir\z -na\= J2 (-irk - na\ + $^(-1)" 

n=-N n=~N n=l 

and by changing n by —n in the first sum we have 

N N N 

\z — na\ = \z + na\ + \z — na\, 

n=—N n=0 n=l 

then the integral given in equation (fT9l) reduces to 



z — na\ 



(20) 



(21) 



Na 



exp 



TV 



N 



a 



\z — na\ 



,n=0 



n=l 



exp 



Na 



N 



N 



^(J2{-ir\z+na\ + Y,i-ir 



— i L2 — "/^a 



dz. 



vra=0 



n=l 



(22) 



The result of the integral given above is (the evaluation of this integral is given in appendix lAl) 



2\A\ 



a 



exp 



:i + 2Ar)aa(-l)^ 



N-l 



sinh 



k=0 



aa 



2eoVo 



eoVo 



a{-l) 



N 



exp 



eoVo 



(23) 



Na 



The result given in equation (|23] ) will diverge unless a — > (— for a > 0. This restriction 
means that we need to have positive infinite charge sheets on both ends of the configuration in 
order for the wave function to go to zero when z — )■ ±oo. Therefore, the normalization constant is 
given by 



A 



-1 



heoVo 


aNa' 


/ exp 




/ ^ 



, 4eoK) 

H exp 

a 



;i + 2N)aa 
2eoVo 



\N 



sinh 



aa 



2eoVo 



where we have used the fact that 

N-l 



-1)^ J2 sinh 



fc=0 



-1 



^k+N 



aa 
2^0 



sinh 



aa 
2^0 



forN> 1. 
(24) 

(25) 



Therefore, if we make the appropiate substitutions we have the bound state energy and wave 
function for a quantum particle of mass m moving through the potential energy function of a 
one-dimensional ionic crystal^, i.e. 



TV 



(26) 



n=-N 



which is given by 



N 

n 



am(-l)"+^ 

n2 



z — na\ 



ri=-N 



(27) 



fi2 r 
— exp - 

ma ^ L 



2maNal 



1 + — exp 

J ma ^ 



{l+2N)ama 
f? 



sinh [^] 



for > 1 and with energy E = — This result shows that it is possible to trap a single electron 
with a one dimensional ionic crystal modeled by repulsive and attractive delta functions for the 
negative and positive ions, respectively. The wave function ip(z) is plotted in Fig. ([U) for different 
values of A^. With the wave function given by equation (1271 ) we can calculate the expectation value 
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FIG. 1. Plots showing the wave function for the following values: (a) = 1, (b) = 2, (c) = 3 and 
(d) = 4. We used atomic units, i.e. h = m = I, and aa = 1. 

of the potential energy, which is given by 



N 



/oo 
exp 

n=-iv 



2am 



TV 



5^(-l)^-+^|z-ja| 



j=-N 



6{z — na)dz. (28) 



solving the integral given above (the evaluation of this integral is given in appendix |B]), we find 
the expectation value of the potential energy function 



ma 



1 + 2A^e-™"''/^' smh[maa/h'^ 



i-(-i)^ 



sinh.[maa / K^] 



Tl? \ ]^ _j_ '2(}—maa/h'^ 

and the average kinetic energy which is given by 

ma^ / 1 1 + 2Ne-""^^/^^ smh\maa/h^] 



(29) 



(T) 



^2 2 



smh[maa / ffi] 



(30) 



Now we will consider the case when the volume charge density is given by 

p = a {5{z — a) + 5{z + a)) , 



the electric field for this charge configuration is given by 



(31) 



and the electrostatic potential will be 



V{z) 



—Cz if z > a 

if —a < z < a 

——Cz if z < —a, 



——z if 2; > a 

-— if -a < z < a 

—z if z < —a, 



(32) 



(33) 



therefore, the dual quantum problem will have a ground state energy given hy E = —2'ma^/ti? 
with the following wave function 

^^-2mc^z/n^ if ;Z > a 

ij{z) = I Ae-2™W?i' if-a<z<a (34) 

J^^2maz/h^ if ^ < -a, 



for the following quantum potential 



U{z) = -a {5{z - a) + 5{z + a)) - 



2ma' 



■9{\z\ < a), 



(35) 



where 9{\z\ < a) is one for —a < z < a and zero otherwise. Notice that the origin of the last term 
in Eq. (l35l) is due to the difference in electromagnetic energy between the asymptotic region and 
all space points. 
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IV. CONCLUSIONS 



We have shown a relation between Poisson and Schodinger equation which enable us to find the 
ground state energy and eigenfunction of a dual quantum system using electrostatic solutions. We 
have obtained the general solution of the wave function for a quantum particle moving through 
a potential energy function that describes a one-dimensional ionic crystal of finite length using 
electrostatic solutions. The result shows that forces exerted by regularly spaced positively and 
negatively ions described by attractive and repulsive Dirac delta functions admit only one bound 
state, i.e. it is possible to trap a single electron in a ionic crystal. The method can be used for non 
periodic potentials which can describe one dimensional ionic quasicrystals. 
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Appendix A: Derivation of the normalization constant 



We need to evaluate the integral 



Na 



Jo J Na 



where 



a 



N 



N 



Y^{-ir\z + na\+Y.^-iy 



\z — na\ 



vn=0 



n=l 



Using the fact that 



(Al) 



(A2) 



a 



N 



n=k+l 



for ka < z < {k + l)a, (A3) 



then we can rewrite the first integral in equation (|AH) in the form 



Na 



^-1 f(k+l)a 



k=Q 



exp 



ka 



a 



eoV'o 



N 



Afz^ (-1)"(2 



na) 



Integrating equation (|A4I) we get 



a — ' 

A:=0 



N 



n=k+l 



aa(-l)''k' 



dz. (A4) 



exp 



n=k+l 



1 — exp 



o"a(— 1) 



(A5) 
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Evaluating the sum over n and regrouping terms in k we end up with 



f(z)j 2eoVo 
e-' ^ 'dz = exp 





aa(-l)^(l + 2iV) 
2^0 



sinh 



fc=0 



aa(—l] 



2eoVo 



In orde to evaluate the second integral in equation (|AH) . we use the fact that 



N 



+ 5^(-ir(2. 



n=l 



for Na < z < oo. 



therefore, evaluating the sum over n we end up with 
integrating equation (IA8I) we have 

eoVo 



oo poo 

e^^^^dz = / exp 

Na J Na 



-z 



Na 



ai-l 



\N 



exp 



eoVo 



dz. 



Na 



(A6) 



(A7) 



(A8) 



(A9) 



This integral given in equation (IA9I) clearly diverges unless we choose a = (— l)^cr, with a > 0. 



Appendix B: Derivation of the potential energy expectation value 



We want to calculate the expectation value of the potential energy, which is given by 



(U) = -a 



N poo 

5^(-ir^|Ar/ exp 

_ AT J —OO 



n=-N 



2am 



N 



j=-N 



S{z — na)dz, (Bl) 



where \A\is the normalization constant. Integrating equation (IB II) we have 



N 



{U) = -a\A\^ (-l)"^'^exp 



n=-N 

by using the following identity 

N n-1 



2maa 



N 



n-j\ 



j=-N 



(B2) 



j=-N 



j=~N 



N 
j=n+l 



n] 



(B3) 



and if we let n — j = k and j —n = I'm the first and second sum, respectively. These substitutions 
change the limits of the sums, and we get 



N 



N-n 



em: 



n-j\ 



(B4) 



j=-N 



k=n+l 



1=1 
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Evaluating the sums over k and / in equation (IB4I) we get 



N 



j=-N 



(B5) 



Substituting equation (IB5|) into equation (IB2I) we have 



(U) = -a|y4|^exp 
Using the fact that 



maa{l + 2N) 



N 



n+N 



exp 



n=-N 



(B6) 



N 



J](-ir+^exp 



n=-N 



maa 



2N 



-1) 



n+N 



exp 



n=0 



maa 



gmaa/?.2^2iVsinh[maa/n2 
(B7) 



and substituting the normalization constant, we end up with 

ma^exp[-^{l + 2N)] fe'^Wft' + 2N smh[maa/h^' 



which can be simplified into 



{l+2N)ama 




ri-(-l)iv] 


fi2 




2 



sinh [^] 



(B8) 



ma 



1 + 2A^e-"'""/^' smh[maa/ff 



^ _j_ 2^Q—m,aa/h'^ 



sinh [ma a/ /i^ 



(B9) 



By using conservation of energy we can obtain the expectation value for the kinetic energy of the 
system. 
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